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Q I Abstract 
(N ■ 

■ We study the effect of additive noise on integro-differential neural field equations. 
Q ! In particular, we analyze an Amari-type model driven by a Q- Wiener process and fo- 

pLn I cus on noise-induced transitions and escape. We argue that proving a sharp Kramers' 

^SJ ■ law for neural fields poses substanial difficulties but that one may transfer techniques 

CN , from stochastic partial differential equations to establish a large deviation principle 
(LDP). Then we demonstrate that an efficient finite-dimensional approximation of the 

■ stochastic neural field equation can be achieved using a Galerkin method and that the 
Oh . resulting finite-dimensional rate function for the LDP can have a multi-scale structure 

' in certain cases. These results form the starting point for an efficient practical compu- 

' tation of the LDP. Our approach also provides the technical basis for further rigorous 
study of noise- induced transitions in neural fields based on Galerkin approximations. 

'. 1 Introduction 

\^ • Starting from the classical works of Wilson/Cowan [56] and Amari [1] there has been consid- 

■ erable interest in the analysis of spatio-temporal dynamics of meso-scale models of neural ac- 
. tivity. Continuum models for neural fields often take the form of nonlinear integro-differential 

equations where the integral term can be viewed as a nonlocal interaction term; see [3l] for 
a derivation of neural field models. Stationary states, travelling waves and pattern forma- 
tion for neural fields have been studied extensively see e.g. [201 128] or the recent review by 
Bressloff [12] and references therein. 
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^ . In this paper we are going to study a stochastic neural field model. There are several 

motivations for our approach. In general, it is well-known that intra- and inter- neuron 
[26] dynamics are subject to fluctuations. Many meso- or macro-scale continuum models 
have stochastic perturbations due to finite-size effects [35| [55]. Therefore, there is certainly 
a genuine need to develop new techniques to analyze random neural systems [16]. For 
stochastic neural fields there is also the direct motivation to understand the relation between 
noise and short-term working memory [48] . 

However, there is relatively little - but fairly recent - work on stochastic neural fields 
which we briefly review here. Brackley and Turner [TT] study a neural field with a gain 
function which has a random firing threshold. Fluctuating gain functions are also considered 
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by Coombes et al. [21]. Bressloff and Webber [15] analyze a stochastic neural field equation 
with multiplicative noise while Bressloff and Wilkinson [TB] study the influence for extrinsic 
noise on neural fields. In all these works the focus is on the statistics of travelling waves 
such as front diffusion and the effects of noise on the wave speed. Hutt et al. [38] study 
the influence of external fluctuations on Turing bifurcation in neural fields. Kilpatrick and 
Ermentrout [IQ] are interested in stationary bump solutions. They observe numerically a 
noise-induced passage to extinction as well as noise-induced switching of bump solutions and 
conjecture that "a Kramers' escape rate calculation" [101 p.l6] could be applied to stochastic 
neural fields but they do not carry out this calculation. In fact, to the best of our knowledge 
there seems to be no general Kramers' law or large deviation principle (LDP) calculation 
available for continuum neural field models although large deviations have been of recent 
interest in neuroscience applications [HI [31]. It is one of the main goals of this paper to 
provide the basic steps towards a general theory. 

Although Kramers' law [5] and LDPs [32] [25] are well understood for finite-dimensional 
stochastic differential equations (SDEs), the work for infinite-dimensional evolution equations 
is much more recent. In particular, it has been shown very recently that one may extend 
Kramers' law to certain stochastic partial differential equations (SPDEs) [H [HI [Z] driven 
by space-time white noise. The work of Berglund and Gentz [7] provides a quite general 
strategy how to 'lift' a finite-dimensional Kramers' law to the SPDE setting using a Galerkin 
approximation due to Blomker and Jentzen [8]. Since the transfer of PDE techniques to 
neural fields has been very successful, either directly [17] or indirectly [221 112], one may 
conjecture that the same strategy also works for SPDEs and stochastic neural fields. 

In this paper we consider a rate-based (or Amari) neural field model 



driven by a Q- Wiener process W for a trace-class operator Q and nonlinear gain function /; 
the technical details and definitions are provided in Section [2l The motivation for ([T|) is given 
in Section [3] where a formal calculation shows that a space-time white noise perturbation of 
the gain function in a deterministic neural field leads to ([T]). 

Remark 1.1. To avoid confusion we alert readers familiar with neural fields that the non- 
linear gain function / in ([1]) is sometimes also called a 'rate function'. However, we reserve 
'rate function' for a functional, to be denoted later by /, arising in the context of an LDP 
as this convention is standard in the context of LDPs. 

Observe that ([1]) is a relatively general formulation of a nonlocal neural field. Hence we 
expect that the techniques developed in this paper carry over to much wider classes of neural 
fields beyond ([T|) such as activity-based models. 

In Section [H we briefly describe important features of the deterministic dynamics for ([T]) 
where e = 0. In particular, we collect several examples from the literature where the classical 
Kramers' stability configuration of bistable stationary states separated by an unstable state 
occurs for Amari- type neural fields. In Section [5] we introduce the notation for Kramers' law 
and LDPs and state the main theorem on finite-dimensional rate functions. In Section [6] we 
argue that a direct approach to Kramers' law via 'lifting' for ([T]) is likely to fail since the 
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Amari model only has a hidden energy-type structure but is not an exact gradient system. 
This raises doubt whether a Kramers' escape rate calculation can actually be carried out 
i.e. whether one may express the pre-factor of the mean first-passage in the bistable case 
explicitly. Based on these considerations we restrict ourselves to just derive an LDP. In 
Section [7] the LDP is established by a direct transfer of a result known for SPDEs. The 
disadvantage of this approach is that the resulting rate function is difficult to calculate, 
analytically or numerically, in practice. Therefore, we establish in Section |H] the convergence 
of a suitable Galerkin approximation for ([T]). Using this approximation one may apply results 
about the LDP for SDEs which we carry out in Section [91 In this context, we also notice that 
the trace-class noise can induce a multi-scale structure of the rate function in certain cases. 
The last two observations lead to a tractable finite-dimensional approximation of the LDP 
and hence also an associated finite-dimensional approximation for first-exit time problems. 
We conclude the paper in Section [TU] with implications of our work and remarks about future 
problems. 



2 Amari-Type Models 

In this study we consider stochastic neural field models with additive noise of the form 



dUAx) 



aUt{x)+ / w{x,y)f{Ut{y))dy 



B 



dt + edWt{x) (2) 



for X G S C R*^, a small parameter e > 0, and t > 0, where B is bounded and closed. In 
([2]) the solution U models the averaged electrical potential generated by neurons at location 
X in an area of the brain B. Neural field equations of the form ([2]) are called Amari-type 
equations or a rate-based neural field models. The equation is driven by an adapted space- 
time stochastic process Wt{x) on a filtered probability space (fi, J-", (J-t)t>o, IP)- The precise 
definition of the process W will be given below. 

The parameter a > is the decay rate for the potential, w : B x B ^ is a. kernel that 
models the connectivity of neurons at location x to neurons at location y. Positive values 
of w model excitatory connections and negative values model inhibitory connections. The 
gain function / : M — M_|_ relates the potential of neurons to inputs into other neurons. 
Typically, the gain functions are chosen sigmoidal and often used are, for example, (up to 
afiine transformations of the argument) f{u) = (1 -|- e~")~^ or f{u) = (tanh(u) + l)/2. 
These examples of gain functions are bounded, infinitely often differentiable with bounded 
derivatives. However, throughout the paper we only make the standing assumption that 

(HI) the gain function / is globally Lipschitz continuous on M. 

We may transfer equation ([2]) into the Hilbert space setting of infinite-dimensional stochas- 
tic evolution equations [23], [50] for the Hilbert space L'^{B). Subsequently brackets (■, ■) 
always denote the inner product on this Hilbert space. Moreover, we introduce the fol- 
lowing notation: Firstly, F denotes the nonlinear Nemytzkii-operator defined from /, i.e., 
F[g){x) = f{g{x)) for any function g G LP'iB) which due to (HI) is a Lipschitz continuous 
operator on L'^{B). Often, also spatially continuous solutions to ([2]) are of interest and thus 
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we note that the Nemytzkii-operator also preserves its Lipschitz continuity on the Banach 
space C{B) with its norm \\g\\o = sup^^^\g{x)\. Secondly, the linear operator K is the 
integral operator defined by the kernel w 

Kg{x)= [ w{x,y)g{y)dx M g E L\B). (3) 

Throughout the paper we assume that 

(H2) the kernel w is such that is a compact, self-adjoint operator on L'^iB). 

We note that an integral operator is self-adjoint if and only if the kernel is symmetric, 
i.e., w{x,y) = w{y,x) for all x,y E B. A sufficient condition for the compactness of K 
is, e.g., ||w;||2,2(gxB) < c>o in which case the operator is called a Hilbert-Schmidt operator. 
If furthermore B is bounded, then the continuity of the kernel w on B x B implies the 
compactness of K considered an integral operator on C{B). 

Then we re-write equation ([2]) as an Hilbert valued stochastic evolution equation 

dUt = [-aUt + KFiUt)] dt + edWt , (4) 

where W is an L^(;B)-valued stochastic process. Interpreting the original equation in this 
form we now give a definition of the noise process assuming that 

(H3) is a Q- Wiener process on L^(S), where the covariance operator Q is a non-negative, 
symmetric trace class operator on L'^{B). 

For a detailed explanation of a Hilbert valued Q- Wiener process and its covariance operator 
we refer to, e.g., [23l [50]. As the operator Q is non-negative, symmetric and of trace class 
there exists an orthonormal basis of L'^{B) consisting of eigenfunctions Vi and corresponding 
non-negative real eigenvalues which satisfy Xli^i Af < oo . It then holds that the Q- 
Wiener process W satisfies 

oo 
i=l 

where /3* are a sequence of independent scalar Wiener processes (cf. [50i Prop 2.1.10]). The 
series ([5]) converges in the mean square on C{[0,T], L'^{B)). Furthermore, a straightforward 
adaptation of the proof of [50l Prop 2.1.10] shows that convergence in the mean-square also 
holds in the space C{[0,T],C{B)) for every T > if G C{B) for all i (corresponding to 
non-zero eigenvalues) and sup^gg Ei^i ^I'^ii^^l < 

The existence and uniqueness of mild solutions to (jlj) with trace class noise for given 
initial condition Uq G L'^iB) is guaranteed under the Lipschitz condition on /, cf. [51], and 
we can write the solution in its mild form 

= e-"*?7o + f e-"(*-') KFiU^) ds + [ e'^^*"") dWs. (6) 
Jo Jo 

The solution possesses a modification in C([0, T], L'^{B)) and from now on we always identify 
the solution with its continuous modification. It is worthwhile to note that for cylindrical 
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Wiener processes - and thus in particular space-time white noise - there does not exists a 
solution to (jl]). This contrasts to other well-studied infinite-dimensional stochastic evolution 
equations, e.g., the stochastic heat equation. Due to the representation of the solution 
([6]) it follows that a solution can only be as spatially regular as the stochastic convolution 
Jq e~"^*~*^ dWg. In the present case the semigroup generated by the linear operator is not 
smoothing in contrast to, e.g., the semigroup generated by the Laplacian in the heat equation. 
Thus the stochastic convolution is only as smooth as the noise which for space-time white 
noise is not even a well-defined function. Particularly, for cylindrical Wiener noise the series 
representation of the stochastic convolution (cf. see equation ([8]) below) does not converge 
in a suitable probabilistic sense. 

We next aim to strengthen the spatial regularity of the solution which will be required 
later on. According to pSf Thm. 7.10] the solution ([6]) is a continuos process taking values 
in the Banach space C{B) if the initial condition satisfies Uq G C{B), the linear part in 
the drift of (jl]) generates a strongly continuous semigroup on C{B), the non-linear term 
KF is globally Lipschitz continuous on C{B) and, finally, if the stochastic convolution is 
a continuous process taking values in C{B). It is easily seen, that the first conditions are 
satisfied and sufficient conditions for the latter property are given in the following lemma. 

Lemma 2.1. Assume that the orthonormal basis functions Vi are Lipschitz continuous with 
Lipschitz constants Li such that 



sup 



< OO 



i=l 



1=1 

for a p E (0, 1). Then the process 

0{x,t):= / e-'^^'-'UWtix) = J2>'i e-^'^'-'Up, 
Jo io 



Y.\^I^'<^ (7) 



possesses a modification with 'y-Holder continuous paths in IR+ x B for all 7 G (0,p/2). 

Proof. We prove the lemma applying the Kolmogorov-Centsov Theorem (cf. [2^ Thm. 3.3]). 
Throughout the proof C is some finite, constant which may change from line to line but 
is independent of x,y E B and t,s > 0. We start showing that the process O is Holder 
continuous in the mean square in each direction. As Vi are assumed continuous these are 
pointwise uniquely given and each 0{t,x) is for fixed x E B and t > a Gaussian random 
variable due to Yli^i '^'i'^ii^)^ < Hence, for all < s < t and all x,y E B we obtain 

E|0(x, t) - 0{y, t)\' = J2 A- / e-2"(*-^) ds \v,{x) - v,{y)\' < C E ^'^'"'^ " ^1' 

1=1 -^^ 1=1 

as \vi{x) — Vi{y)\'^ < Lf''\x — y\''\x — y\'' for every p G [0, 1] following [231 Sect. 5.5.1] and 
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\x — < diam(i3)''. Next, for the temporal regularity we obtain 
E\0{x,t) -0{x,s)\'^ 

CO „± oo 

-a(t-r) _ g-o(s-r) |2 



^A^^;,(x)2 / e-2"(*-^)dr + ^A2^;,(x)2 / 
i=i i=i 

°° /p-aO _ „-a{t-s) \ / fp-aO _ „-a{t-s)\2 _ ( -at _ -as\2\ 

<E^M^r {- ^ )+cj::qv.(xf ((^ ^ ^) 



As the exponential function is on the negative half axis Holder continuous for every p G [0,1] 
it holds 

oo 

E\0{x,t) -0{x,s)\^ < C^A,2^;,(x)2|^-s|^ 

i=l 

Thus, overall Jensen's inequality yields E|0(x, t) — 0{y, s)p < C{\x — y\'^+ \t — s\'^)^^'^. Since 
the difference 0{x,t) — 0{y,s) is centered Gaussian it further holds that 

E\0{x,t) -0{y,s)\^"' < C{\t- + \x-y\^r''/^ VmGN. 

Now, the Kolmogorov-Centsov Theorem implies the statement of the lemma. □ 

We present an example to illustrate the type of noise we are generally interested in. 
Further motivation is provided in Section |3l 

Example 2.2. Consider the neural field equation on a d- dimensional cube B = [0,2'itY with 
noise based on trigonometric basis functions of L^{[0,27rY). This type of noise is almost 
ubiquitous in applications as for the stochastic heat equations the basis functions can be cho- 
sen such that the usual (Dirichlet, Neumann or periodic) boundary conditions are preserved. 
For the example of noise preserving homogeneous Neumann boundary conditions the basis 
functions are 

d 

Vi{,x) = Wei^{xk), (9) 

k=l 

where x = (xi, . . . , Xd), i = {ii, . . . ,ik) is a multi-index in N'^ and the functions are given 
by 

1 

«fc = 0, 



V27r 

—= cos{ikXk/2) ik>l- 



The functions Vi are for all i E N'^ pointwise bounded by tt"*^/^ and Lipschitz continuous 
with Lipschitz constants given by Li = n"'^^'^ \i\ (cf. /2l Lemma 5.3]). Next, we construct 
a trace class Wiener process from these basis functions. A particular important example of 
spatio-temporal noise is smooth noise with exponentially decaying spatial correlation length 

1 / Ti \x — yp\ 

EWt{x)Ws{y) = minjt, s} -— — -7 exp ( — — — -- — ) + correction on the boundary (10) 

(24)'' V 4 c, ^ 
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for a parameter ^ > modelling the spatial correlation length. Note that for ^ — > this 
noise process approximates space-time white noise. Following we can calculate under 
the assumption that ^ <^2tt the coefficients Af such that the Q- Wiener process possesses 
the correlation function ffTOj) and obtain 



A 



exp 



4tt 



Now, it is easy to see that for this choice of eigenvalues the noise is of trace class and moreover 
the addtional conditions of Lemma \2.1\ are satisfied: As the functions Vt are bounded we obtain 



sup 



< TT + TT 



Yl Yl 

N=l ie{0,...,N}'i\{0,...,N-l}d 



2 L-|2 



Ait 



< 71 +71 



E 

N=0 



exp 



< OO 

and the second condition of ([7]) is satisfied as 

oo 

EAji?" < T^'^E E 

i=N'* N=l i(^{0,...,N}''\{0,...,N-l}<' 



exp 



An 



|2p 



'^E 



exp 



N=0 



Att 



N-l 



< oo , 



3 Gain Function Perturbation 

Another motivation for the considered additive noise neural field equations stems from a 
(formal) perturbation of the gain function / with space-time white noise. Let VV denote 
space-time white noise and consider the randomly perturbed Amari equation 

dtU{x, t) = -aU{x, t) + w{x, y) (^f{U{y, t)) + e Mv, t)) dy . (12) 

Recall that, by assumption (H2), the integral operator K defined by the kernel w is a self- 
adjoint compact operator. Thus the Spectral Theorem implies that K possess only real 
eigenvalues Aj, i G N, and the corresponding eigenf unctions Vi form a basis of L'^{B). If 
additionally we assume that 

(H4) i^' is a Hilbert-Schmidt operator on L?{B), that is, ||w||L2(gx5) < oo, 

then the eigenvalues satisfy Yl^i -^f < Hence, K possesses the series representation 

oo 
i=l 
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which yields for the perturbed equation (!T2|) the representation 

dtU{x, t) = -aU{x, ^) + £ f(U{y, t))v,{y) dy + e (>V(t, ■), tk)^ tk{x) . 

Next, note that the random variables PI = (>V(-, t), Wj) form a sequence of independent scalar 
white noise processes in time. Therefore, the perturbed equation becomes 

„ oo 

dtU{x,t) = -aU{x,t) + / w{x,y)f{U{y,t))dy + ey2^i^>ii^)- 

i=i 

Re-writing this equation in the usual notation of stochastic differential equations we obtain 
dUt{x) = -aUtix) + I w{x,y)f{Ut{y))dy + edWt{x), (13) 



where 

oo 

Wt{x) = Y,^^PK^) 
i=l 

is a trace-class Wiener process on the Hilbert space L'^iB). Note, when comparing to ([5]) 
here the coefficients Aj may be negative, however, as also — /3* is also a Wiener process this 
slight inconsistency can be neglected. 

We next want to discuss spatial continuity of the solution to this equation with its 
particular noise structure. It is clear that this should translate into smoothing conditions of 
the kernel w. Due to Lemma [2.11 it is sufficient to establish conditions ([7]): First, it holds 
that 

oo oo ^ „ \ 2 oo 

VA-Wi(x)^ = V( / w{x,y)vi{y)dy\ =y^{w{x,-),Vif = \\w{x,-)\\l2 
i=l i=i / i=i 

due to Parseval's identity. Hence, the first condition of ([7]) becomes 

sup ||w(x, ■) 11^2 < oo. (14) 

Next, the basis functions are continuous if the kernel w{x, y) is continuous in x and as the 
minimal Lipschitz constant is given by the supremum on the derivatives we obtain 



sup 



-^V^ / w{x,y)vi{y)dy 
Jb 



< -^sup||V^w(x, ■)||i2 



xGB 

due to the Cauchy-Schwarz inequality. Therefore the second condition in (j7]) is satisfied if 



oo 



sup ||Va;M;(a;, ■)||i2 < oo and ^ < oo (15) 



XGB .^^ 



for a p G (0, 1). The condition ( !T^ and the first part of (fT5|) are easily checked but for the 
second part of (fT5l) usually theoretical results on the speed of decay of the eigenvalues have 
to be obtained. We note that ( ITSi) is certainly satisfied if is a trace class operator with 
P = 1/2. 
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4 Deterministic Dynamics 

The classical deterministic Amari model, obtained for e = in ([2]), is 

dtU{x, t) = -aU{x, t)+ [ w{x, y)f{U{y, t))dy. (16) 



B 



where B C W^. Suppose there exists a stationary solution U* = U*{x) of (fT6|) . To determine 
the stability of U* consider U{x,t) = U*{x) + %lj{x,t). Substituting into ( !T6|) and Taylor- 
expanding around U* yields the linearized problem 



dMx.t) = -aijix,t) + / wix,y){Df){U*{y))^lj{y,t)dy. (17) 

Jb 

Hence, the standard ansatz ip^x^t) = ?/;o(x)e^* leads to the eigenvalue problem 

{fi + a)M^)= [ w{x,y){Df){U*{y))My)dy:={CtlJo){x) or £^0 = P^o- (18) 



B 



-.7] 



The linear stability condition /i < is equivalent to rj < a where p G spec(£). The stability 
analysis can be reduced to the understanding of the operator C. However, this is a highly 
non-trivial problem as the result depends upon B, U*{x), w{x,y) and f{u). 

An LDP and Kramers' law are of particular interest in the case of bistability. Therefore, 
we point out that there are many situations where (fT6|) does have three stationary solutions: 
U^{x) which are stable and Uq{x) which is unstable. The following three examples make 
this claim more precise. 

Example 4.1. The first example is presented by Ermentrout and McLeod f28^ . Let B = M., 
w{x,y) = w{\x — y\), a = 1 and assume that < U{x,t) < 1. Furthermore, suppose that 
f E Ci([0,l],M) with /' > and 

f{U) ■.= -U + f{U) 

has precisely three zeros U = 0, a, 1 with < a < 1. The additional conditions /'(O) < 1 and 
f'{l) < 1 guarantee stability of the stationary solutions U = and U = 1. As an even more 
explicit assumption 12^ p. 463] one may consider a Dirac 6 -distribution for w in ( IT6|) which 
yields 

dtU{x,t) = -U{x,t) + F{U{x,t)). (19) 

The stability condition for (I19p is F'{U*) < 1 and F'{a) > 1 yields an unstable stationary 
solution between the two stable states. 



Example 4.2. An even more concrete example is given by Guo and Chow 13^ \3^ . They 
assume i3 = M, w{x,y) = w{x — y) , a = 1 and fix two functions 

f{u) = [b{u - Ub) + l]H{u - Ub), w{x) = Ae-'^l^l - e'^' 

where H{-) is the Heaviside function and h, a, A and Ub are parameters. Depending on 
parameter values one may obtain three constant stationary solutions exhibiting bistability as 



expected from Example 4-F However, there are also parameter values so that three stationary 



pulses exhibiting bistability exist. 
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Examples I4.HI^I2] are typical for many similar cases with x G M or x G M^. Many results 
on existence and stability of stationary solutions are available see e.g. [H HTJ |131 US] and 
references therein. 

Example 4.3. As a higher- dimensional example one may consider the work by Jin, Liang 
and Peng fS^ who assume that w{x, y) = w{x — y), a = 1, B = and 

Zoo = / w{x) dx < oo, K,Woo > 1- 

where k is the Lipschitz constant of f E C^(]R'',M). Furthermore, suppose f is uniformly 
continuous and 

f'iU)Zoo < 1 /orf/G (-oo,t/i)U(t/2,oo), 
f'{U)Zo^ = 1 forUe{Ui,U2}, 
f'{U)Zo, > 1 forUeiUuU2), 

for Ui < < 1/2- Then fSR Prop. 11] the conditions 

-Ui + /(f/i)Zoo < and -U2 + f{U2)Zoo > 

yield three stationary solutions U^, f/^ and Uq. The solutions f/^ are stable and satisfy 
UI <0 and [/| > 0. The solution Uq is unstable. 

Although we only focus on stationary solutions, it is important to remark that the tech- 
niques developed here could - in principle - also be applied to travelling waves U{x, t) = 
U{x — st) for s > 0. The existence and stability of travelling waves for (fT6l) has been investi- 
gated for many different situations; see e.g. [121 1211 [ISl [22] and references therein. However, 
it seems reasonable to restrict ourselves here to the stationary case as even for this simpler 
case an LDP and Kramers' law are not yet well understood. 



5 Large Deviations and Kramers' Law 

Here we briefly introduce the background and notation for LDPs and Kramers' law needed 
through the remaining part of the paper; see [321 125] for more details. Consider a topological 
space X with Borel a- algebra Bx- A mapping / : A" — )■ [0, 00] is called a good rate function 
if it is lower semicontinuous and the level set {h : I{h) < a} is compact for each a G [0, 00). 
Sometimes the term action functional is used instead of rate function. Consider a family 
{/i''} of probability measures on {X, Bx)- The measures {if} satisfy an LDP with good rate 
function I if 

- inf / < lim inf In /i' (F) < lim sup In /i' (F) < - inf / (20) 

r° e^o r 

holds for any measurable set of paths T d X\ often infima over the interior F° and closure 
F coincide so that lim inf and lim sup coincide at a common limit. One of the most classical 
cases is the application of fl2Q|) to finite-dimensional SDEs 

dut = g{ut)dt + eG{ut)d(3t (21) 
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where Ut G M^, ^ : -> M^, G : ^ M^^'^ and A = • • • , PtY is a vector of k 
independent Brownian motions. The topological space is chosen as a path space 

X := Co([0, T],M^) = {0 G C([0,T],M^) : 0(0) = 0}. 

To state the next result we also need the Sobolev space 

Hi ■= {0 : [0, T] ^ : absolutely continuous, 0' G L^, 0(0) = 0}. 

Furthermore, we are going to assume that the diffusion matrix D{u) := G{uf-G{u) G M^^^ 
is positive definite. 

Theorem 5.1. (JS^,JE3) The SDE ^ satisfies the LDP ^ given by 

- inf / < lim inf In ¥iut G T) < lim sup In ¥iut G T) < - inf /. (22) 

r° e^o e^o r 

for any measurable set of paths T d X with good rate function 

, / i/o'(0;-/(0,))^2)(0*)"H01-/(0t))rft, ^/0e«o + ^ri, 

![(()) = I[o,T][(t)) = { , (23) 

[ +00 otherwise. 

Probably the most important application of the LDP (l22l) is the so-called first-exit prob- 
lem. Suppose that ut starts near a stable equilibrium u* E V C where P is a bounded 
domain. Define the first-exit time 

TV := inf{t >0:ut^V}. (24) 

To formalize the application of the LDP define the mapping 

Z{u, v; s) := inf {J(0) : G C([0, s], M^), 0o = 0, = v} (25) 

which is the cost for a path starting at u to reach v in time s. By Theorem 15.11 it follows 
that 

lime^lnP(ri, < t\u = Uq) = M{Z{u,v;s) : s E[0,t],v ^ V}. (26) 

For the case of a stable deterministic equilibrium u* G V, where T> is inside the basin of 
attraction of u*, the function 

Z(u*, v) = inf Z(u*, v\ t) 

is called the quasipotential. It is natural to minimize the quasipotential over &D and define 

Z := inf Z(u*, v). 
yedv 

Theorem 5.2. (f ^ .fE ^ ) For all initial conditions u eT) and all 5 > 0, the following two 
limits hold: 

hmP fe(^-^)/^' <Tv< e(^+^)/^Vo = u*) = 1, (27) 
\\me^\nW.[T'D\uQ = u*] = Z. (28) 
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If the SDE (I2T]) has a gradient structure with identity diffusion matrix i.e. 



g{u) = -VV{u)hTV:W'->R, and G{u) = IdeW''' (29) 

then one can show [32] that the quasipotential is given by Z{u*,v) = 2{V{v) — V{u*)). If 
the potential has precisely two local minima and a saddle point u* with N — 1 stable 
directions so that the Hessian V'^V{u*) has eigenvalues 

Pl{ul) < < P2«) < ■■■ <P7V«) 

then one can even refine Theorem 15.21 Suppose Uq = u*_ then the mean first passage time 
to u*^ satisfies 



.M,. > : 11. -u-J.<_ ^ ^J^^^,^n..-n.u.. ,30) 



which is also known as the Kramers' law [5j or Arrhenius-Eyring-Kramers' law [21 [291 112] • 
Note that the key differences to the general LDP (1271) for the first-exit problem are that flHU]) 
yields a precise prefactor for the exponential transition time and uses the explicit form of 
the good rate function for gradient systems. It is interesting to note that a rigorous proof of 
( 130|1 has only been obtained quite recently [9l [TO] . 



6 Gradient Structures in Infinite Dimensions 



The finite-dimensional Kramers' formula (!30|) applies to SDKs ( 12T|) with a gradient-structure 
g{u) = — VI^(m) where V : M is the potential. A generalization of Kramers' law has 

been carried out to the infinite-dimensional case of SPDEs given by form 



dUt = [AU - h'{U)]dt + e dW{x, t) 



(31) 



for U = U{x, t), X G i3 C M, i3 a bounded domain, h E C'^(M, M) for suitably large k eN and 
W{x, t) denotes space-time white noise and either Dirichlet or Neumann boundary conditions 
are used [Hi^l [Z] • A crucial reason why this generalization works is that the SPDE (IM]) has 
a gradient-type structure [30] given by the energy functional 



V[U] : = 



IB . 



-u'ixr+hiuix)) 



dx. 



Indeed, computing the Gateaux derivative yields 



VzV[U] 



-U"{x) + h'{U{x))]z{x) dx. 



(32) 



(33) 



The Gateaux derivative is equal to the Frechet derivative W = DV by a standard continuity 
result [2^ P-47]. Hence (l33ll shows that the stationary solutions of ( pTl) are stationary points 
of the energy functional V. However, for the Amari-model (fT6l) the analogous naive guess 



V[U] :- 



a 



U{x) 




U{y) 



f{r)w{x, y)drdy 



B Jo 



dx 



(34) 



12 



fails direct computation yields 



\\m-{y[U + 5z]-V[U]) 

5--s>0 



all {x)z{x) 



f{U{y))w{x,y)z{y)dy 



dx. 



(35) 



Therefore, setting VzV^[f/] = is not equivalent to the solution of the stationary problem 

-aU{x)+ I w{x,y)f{U{y))dy = 0. 
Jb 

Due to the presence of the different terms z{x) and z{y) in ( l35l) one may guess that the 
modified functional 



V[U] :- 



B L 



2^(^) I f{U{y))f{U{x))w{x,y)dy 



dx 



(36) 



could work. However, another direct computation shows that 



VzV[U] 



[aU{x)z{x)dx] — - 




f{U{x))Df{U{y))w{x,y)z{y)dy dx 




UB JB 

fiU{y))Df{U{x))w{x,y)z{x)dy dx 



'B JB 

= a{U,z)-{KF{U),Df{U)z) 

Hence / and its derivative Df both appear instead of the desired formulation; by a similar 
computation one can show that replacing /(«(■)) in fl36|) by Jq f{r)dr fails as well. Hence 
there does not seem to be a natural gradient formulation. However, one has to consider pos- 
sible coordinate changes. The idea to apply a preliminary transformation has been discussed 
e.g. in [271 P-2] and [SI p.488]. Assume that 

=: g exists and g' ^ 0. (37) 

Define P{x,t) := f{U{x,t)) as the mean action-potential generating rate so that U = g{P)- 
Observe that 

^ -ag{P{x,t)) 



dtP{y,t) 



w{x,y)P{y,t)dy 



(3^ 



g'{P{x,t)) 

For this equation the problem observed in fl36l) should disappear as the integral only contains 
linear terms. One may define an energy-type functional 

.P(x) 



E[P] 



Calculating the derivative yields 



ag{r)dr — - / w{x,y)P{y)P{x)dy 



dx. 



VqE[P] = hm-{E[P + 6Q]-E[P]) 

<5-s-0 



ag{P{x))Q{x)dx [ [ w{x,y){P{y)Q{x) + P{x)Q{y)}dy dx, 
^ Jb Jb 



{ag{P),Q)- 



w{x,y)P{y)dy,Q ] . 
13 



This shows that there is hidden energy-type flow structure in the Amari model for the 
assumptions (1571) so that 

dtP{x,t) = - ] VE[P{x,t)]. (39) 

However, even with this variable transformation there seems to be httle hope to derive a 
precise Kramers' rule for the stochastic Amari model ([2]). The problems are as follows: 

• There is still a space-time dependent nonlinear prefactor l/g'{P{x,t)) in fl39p for the 
deterministic system, so the system is not an exact gradient flow for a potential. 

• Applying the change-of-variable Pt{x) := f{Ut{x)) for the stochastic Amari model ([2]) 
requires an Ito-type formula so that 

dt + eM{Pt{x))dWt{x), 
(40) 

where M(Pt(x)) is now a multiplicative noise term; see [IS] and references therein for 
more details on infinite-dimensional Ito-type formulas. The higher-order term 0{e'^) 
in the drift part of (HOj) is not expected to cause difficulties but a multiplicative noise 
structure definitely excludes the direct application of Kramers' law. 

• Even if we would just assume - without any immediate physical motivation - that the 
noise term in P0|) is purely additive edWt{x) there is a problem to apply Kramers' law 
since we do not have a structure like in fl^Tl) with G{-) = Id as Wt{x) is a Q- Wiener 
process defined in and driving space-time white noise in (jl]) is particularly excluded 
due to the non-existence of a solution. 

Based on these observations an immediate approach to generalize a sharp Kramers' for- 
mula to neural fields seems unlikely. Hence we try to understand an LDP for the stochastic 
Amari-type model ([2]). 

7 Direct Approach to an LDP 

A general direct approach for the derivation of an LDP for infinite-dimensional stochastic 
evolution equations is presented in [23] and further results are have been obtained for certain 
additional classes of SPDEs [TTl [TSl [ISl E2]- The results in [23] are valid for semilinear 
equations with suitable Lipschitz assumptions on the nonlinearity and with solutions taking 
values in C{T>). We state the available results applied to continuous solutions of the Amari 
equation (jl]) assuming that the conditions of Lemma 12.11 are satisfied. 

For the following we assume that there exists an open neighbourhood V G C{B) con- 
taining a stable equilibrium state u* of the deterministic Amari equation flTBl) such that f> 
is contained in the basin of attraction of u*. We are interested in the rate function and the 
first-exit time of the process from V given by 

Tp = inf {t >0:U^V} 



dPAx) 



1 



g'iPtix)) 



-ag{Pt{x) 



w{x,y)Pt{y)dy + 0{e^ 
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if U starts in the deterministic equilibrium state u*. In order to state the quasipotential for 
u we consider the control system 

y = -ay + KF{y) + Q'/^v, y^ = xeC{B) (41) 

for controls v G L^((0, T), L^(,B)) for all T > and denote by y^'"" its unique mild solutior0 
taking values in C([0, T], C{B)) for all T > 0. Then we define 

I{u\z)=mi\^-j^^ \\v{s)\\l,ds: |/"*'^(T) = T > o| , (42) 

where this quasi-potential relates to the minimal energy necessary to move the control system 
f HTj) started at the equilibrium state u* to z. 

Theorem 7.1. (l2^ Thru. 12.18]) It holds that 

lime^ InErr-plf/n = n*l = inf I(u*,z). 

Following further the exposition in [231 Sec. 12] explicit formulae for the rate function I 
are only available in the special case of the drift possessing gradient structure and space-time 
white noise. As we have argued above, this structure is particularly not satisfied for neural 
field equations. Hence, the same observations as presented at the end of the last section 
prevent a further direct analytic approach to the LDP. Therefore, we try to understand the 
LDP problem for a discretized approximate finite-dimensional version of the neural field 
equation. 



8 Galerkin Approximation 

Throughout the section we assume that the assumptions (H1)-(H3) are satisfied. As a dis- 
cretized version of the neural field equation we consider its spectral Galerkin approximations. 
In order to decouple the noise, we define the spectral representation of the solution 



Ut{x) = ^ul Vi{x) . 



(43) 



1=1 



Here the orthonormal basis functions Vi are given by the eigenfunctions of the covariance 
operator of the noise with corresponding eigenvalues A^, see equation To obtain a 



^The existence of such a solution is guaranteed by standard results on deterministic equations Ccf. p3l 
Sect. A. 3]) as long as Q^^^ maps L^{B) continuously into C{B). This is easily established. The unique square 
root Q^/^ of Q is the Hilbcrt-Schmidt operator given by Q^^^g = X^i^i ^ii'^iTg)'Vi for all g G L?{B) and in 
order to show that Q^l'^g £ C{B) it remains to establish that the functions converge uniformly on B. This 
holds as for all x € S 



^A,(«„gK(s)| < |^A^>,(x)2|'^'|^(z;„g)2p'' < (^sup|^A^>,(a;)2|y(^^(u„5) 

i=7V i=N 4=1 '^'^^ i=l i=N 



1/2 



1/2, 



1/2 



Hence the upper bound, which is finite due to (O , is in independent of x and converges to zero for N 
Moreover, we further find that g{t) e L'^{{0,T), L^{B)) implies Q^''^g{t) £ I^^{{Q,T),C{B)). 
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equation for the coefficients we take the inner product of equation (jlj) with the basis 
functions Vi which yields 

{dUt,v,) = [-aiUt,Vi) + {KF{Ut),Vi)]dt + e{dWt,Vi) for t e N. 

After plugging in f l43|) we obtain for the countable Galerkin system 

dui = [-aui + {KFYiul u\, . . .)] dt + e\id^\ for i G N. (44) 

Here the non-linearities coupling all the equations are given by 



j w{x,y) f i^Yy^t Vj{y)^ dy^ dx 



w{x,y)vi{y)dy j dx 

due to the symmetry of the kernel w. If, in addition, we assume that (H4) holds and K and 
Q possess the same eigenfunctions and the eigenvalues are related as discussed in Section [3] 
the non-linearities become 

{KFYiulul...) = ^^jj [fyi vA^)^ v,{x)dx. (45) 

The A^-th Galerkin approximation to U is obtained truncating the spectral representation 
( 143|) and thus given by 



N 



f/f = 5^«;'^^,, (46) 



i=l 



where U-^^ are the solutions to the A^- dimensional Galerkin SDE system 

dv^f = [-aut'' + (J^F)^'^(m^'^, . . . , ''^)] dt + e\id(3l = 1, . . . , A^, (47) 

where the non-linearities KF^'^ are given by 

(J^F)''^(«i'^,...,«^'^) = jj {fy^^"" v,{x)^ (^l^w{x,yMy)dy^dx (48) 

or, in the special case of Section [3l by 

(KF)^'^(«i'^, . . . , w^'^) = \jj{^^ v,{x)^ v,{x)dx , (49) 

respectively. The following theorem establishes the almost sure convergence of the Galerkin 
approximations to the solution of @. Therefore, we may be able to infer properties of the 
behaviour of paths of the solution from the path behaviour of the Galerkin approximations. 
We have deferred the proof of the theorem to the appendix. 
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Theorem 8.1. It holds for all T > that 

lim suptg[o,T] \\Ut - t/f = a.s. 

If, in addition, the series Yl^i^l'^f converges in C{B) and the functions Vi are Lipschitz 
continuous with Lipschitz constants Li such that there exists an p E (0, 1) with Xli^i -^i -^i'' < 
oo (i.e. the conditions of Lemma \2.1\ are satisfied), Uq e C{B) such that hniTv-^oo ||f^o ~ 
P^UqWo = and K is compact on C{B), then it holds for all T > that 

hm supigfo ||f/t — f//^||o = a.s. 



9 Approximating the LDP 

The LDP in Theorem 17.11 is not immediately computable. Here we show that a finite- 
dimensional approximation can be made and what the structure of this approximation en- 
tails. For simplicity, consider the case when the diagonal diffusion matrix D with entries 
2)jj = is positive definite i.e. Aj 7^ for all i G N. Observe that the inverse of D induces 
an inner product on for G N U {00} via 

(a, b)N := (D)^' b = [D~'^^af [1)-^/%] for a, b G . 

where D is understood as the projection onto R^^^ if < 00. We are also going to use the 
notation introduced in Section [8] for the Galerkin approximation i.e. u'^^ denotes the vector 

/ 1,N 2,N N,N^T r- toA^ 

Theorem 15.11 immediately implies the following: 

Proposition 9.1. For the finite- dimensional Galerkin system (1471) the rate function is given 
by 

jN^^^^^ Uo((<f>fy-9'''{<Pf)d<Pfy-9'''{<pf))Ndt, ^f<j>euo + H,, ^^^^ 

foo otherwise, 

where g^'^'i^j^f) = -a<l^''' + (K Ff'^ {<f>]''' , . . . , 

Recall from Section [7] that Theorem 17. II provides a large deviation principle. For the case 
when Q is a positive operator, we may formally re-write the control system ( 14T|) as 

D~'%~{-ay + KF{y))]=v (51) 

so that the rate function for the Amari model can be expressed as 

J(0) = <' Uo l5^~'^'Wt-9m^~'^'Wt-9i<l>t)]dxdt, if0G«o + i^i, ^^2) 
-00 otherwise, 

where g{(j)t) = —a(j)t + KF((j)t) and T)~^^'^u = Therefore, the next result 

just implies that the Galerkin approximation is consistent for the LDP. 
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Proposition 9.2. For each 0j G Mq + Hi we have liiiiTv-i^oo |-^( 



0. 



Proof. Considering tlie finite-dimensional rate function ( 1501) it suffices to notice that 



9' 



N 

E 

i=l 
N 

E 

i=l 



A? 



X2 



by orthonormality of the basis in L'^{B). 



dx 



□ 



Hence we may work with the finite-dimensional Galerkin system and its LDP for compu- 
tational purposes. However, the truncation may still be very large. We are going to show, 
using a formal analysis for a certain case, that there is an intrinsic multi-scale structure of 
the rate function. We assume that we are in the special case considered in Section 12] where 
K and Q have the same eigenfunctions and the corresponding eigenvalues are given by Aj 
and Af, respectively. 

Lemma 9.3. For each N E'N the first part of the rate function (150|1 can he re-written as 



I^((P) = 11 ai - 2a2 + ag dt 



T 



(53) 



where the three terms are given by 



and {KF) 



i,N 



UKF) 



i.N 



N 



{{<pfy + a<pf,KF-'^{<pf)) 



N, 



Proof. For notational simplicity we shall omit in this proof the subscript for the inner product 
(■) ')n = ■) as well as the Galerkin index e.g. (p^^ = (pt as it is understood that we work 
with A^-dimensional vectors. Consider the following general calculation 

{^[-g(^,),^[-g(^,)) = {<p[,<p[)-2{<p[,g{<Pt)) + {g{cl>t),g{cl>t)) 

= {<P[,<P[) + 2a{<P[,<Pt)-2{<P[,KF{<P,)) 

+ {KF{(Pt),KF{^t)) +a^q^t,4>t) -2a{(Pt,KF{(Pt)) 
= {<j)[ + a<j)t, <l)'t + a(j)t) - 2(0; + a<j)t, KF{<f)t)) + K~F{<l)tf K~F{<j)t). 
and observe that the result is independent of A^. □ 

It is important to point out that the LDP from Theorem 15.11 requires the infimum of the 
rate function. From Lemma 19.31 we know that the rate function splits into three terms. The 
three terms are interesting in the asymptotic limit N ^ oo. Suppose 



and 



~ • N 

KF ' ( 



as — oo for some non-negative functions k, 77. Then Lemma [9.31 yields 

af = 0{K{NyX],'), af = 0{K{N)viN)\-^'), = 0{r^{Nf). 
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Lemma 9.4. Suppose there exists a positive constant Kf such that 



snp\fix)\<Kf (54) 



then r]{N) = 1. 

Proof. A direct estimate yields 

N 





<-L 







/(5^0i'%.(x)) 



i=l 



\vj{x)\dx < Kf I \vj{x)\dx. 



Since ||t'j||L2{B) = 1 and L^{B) ^ L^{B) the last integral is uniformly bounded over j G N 
by meas(i3)i/2_ □ 

We remark that several typical functions / discussed in Section [2] such as f{u) = (1 + 
e~")~^ and f{u) = (tanh(u) + l)/2 are globally bounded so that Lemma [9^ does apply to 
many practical cases. In this situation we get that 

af = 0{K{NfX-'), = 0{k{N)\-'), = 0{1). 

We make a case distinction between the different relative asymptotics of k,{N) and Xn- 

• If K,[N) ^ Aat or K,{N) ~ Ajv as — oo then we can conclude that k{N) — i.e. 

{(f)^'^y + a(f)^'^ ^0 as iV ^ (55) 

since for trace-class noise we know that Aat — > 0. If we formally require that {<p^'^)' + 
acf)^'^ = for sufficiently large then the higher-order Galerkin modes decays expo- 
nentially in time 

iN,N iN.N -at 

• If k{N) ^ Aat as A^ — )■ oo then ai ^ — 2a2 + as and the first term dominates the 
asymptotics. But > for all A^ so that the rate function only has a finite infimum 
if af^ — )■ as A^ — )■ oo. This implies again that ( 155|1 holds for the case of a finite 
infimum. 

Hence we get in many reasonable first-exit problems for the Amari model with trace- 
class noise that there is a finite set for n < A^ of 'slow' or 'center-like' directions and an 
infinite set of 'fast' or 'stable' directions for n > N. Although we have made this observation 
from the rate function alone, it is entirely natural considering the structure of the Galerkin 
approximation. Indeed, for the case when the eigenvalues of K and Q are related we may 
write ( 1471) as 

dul'^ = ^— au^'^ + \i[ ■ ■]j dt + eXid(3l (56) 

so that for bounded nonlinearity /, which is represented in the terms [■■■] in ( 156|) . the 
higher-order modes should really just be governed by du]"'^ = —aul'^ dt. 
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Hence, Propositions 19.1119. 21 and the multi-scale nature of the problem induced by the 
trace-class noise suggest a procedure how to approximate the rate function and the associated 
LDP in practice. In particular, we may compute the eigenvalues and eigenf unctions of K 
and Q up to a sufficiently large given order N* . This yields an explicit representation of the 
Galerkin system and the associated rate function. Then one may apply any finite-dimensional 
technique to understand the rate function. One may even find a better truncation order < 
A^* based on the knowledge that the minimizer of the rate function must have components 
that decay (almost) exponentially in time for orders bigger than N . 

10 Outlook 

In this paper we have discussed several steps towards a better understanding of noise-induced 
transitions in continuum neural fields. Although we have provided the main basic elements 
via the LDP and finite-dimensional approximations, there are still several very interesting 
open problems. 

We have demonstrated that a sharp Kramers' rate calculation for neural fields with trace- 
class noise is very challenging as the techniques for white-noise gradient-structure SPDEs 
cannot be applied directly. However, we have seen in Section H] that the deterministic dy- 
namics for neural fields frequently exhibits a classical bistable structure with a saddle-state 
between stable equilibria. This suggests that there should be a Kramers' law with exponen- 
tial scaling in the noise intensity as well as a precisely computable pre-factor. It is interesting 
to ask how this pre-factor depends on the eigenvalues of the trace-class operator Q definin- 
ing the Q- Wiener process. We expect that new technical tools are needed to answer this 
question. 

Another interesting question - which we aim to address in future work - is to use the 
Galerkin approximations in Section |8] and the associated LDP in Section [9] for numerical 
purposes [15]. In particular, recent work for SPDEs [8] suggests that a spectral method can 
also be efficient for stochastic neural fields. Results on numerical continuation and jump 
heights for SDKs can also be immediately transferred to the spectral approximation 
which would allow for parameter-dependent studies of bifurcations and associated noise- 
induced phenomena. 

One may also ask how far the technical assumptions we make in this paper can be 
weakened. It is not clear which parts of the global Lipschitz assumptions may be replaced 
by local assumptions or removed altogether. Similar remarks apply to the multiscale nature 
of the problem induced by the decay of the eigenvalues of Q. How far this observation can 
be exploited to derive more efficient analytical as well as numerical techniques remains to be 
investigated. 

On a more abstract level it would certainly be desirable to extend our basic framework 
to other topics that have been considered already for deterministic neural fields. A gen- 
eralization to activity based models with nonlinearity f{j^w{x^y)u{y)dy) seems possible. 
Furthermore, it may be highly desirable to go beyond stationary solutions and investigate 
noise-induced switching and transitions for travelling waves and patterns. 
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A Convergence of the Galerkin Approximation 

Proof of Theorem \8.1[ We fix a T > 0. Throughout the proof an unspecified norm || • || or 
operator norm ||| ■ |||, respectively, are either for the Hilbert space L'^{B) or the Banach space 
C{B) and estimates using the unspecified notation are valid in both cases. Furthermore, 
C > denotes an arbitrary deterministic constant which may change from line to line but 
depend only on T. We begin the proof obtaining an a-priori growth bound on the solution 
of the Amari equation (jl]). Using the linear growth condition on F implied by its Lipschitz 
continuity we obtain the estimate 

\m\ < e-"*||f/ol| +C re-"(*-^)(l+ ||f/J|)ds+ IIOJ . 



Due to Gronwall's inequality there exists a deterministic constant C such that it holds almost 
surely 

sup \\Ut\\ < C(l + \\Uo\\ + sup \\0t\\) e^^ a.s. (57) 

Note that O is an Ornstein-Uhlenbeck process and it thus holds sup^g^j-] || 0(11^2 < oo 
almost surely and under the assumptions of Lemma |2TT] in addition sup^gjo yj ||0j||o < oo 
almost surely. 

Let denote the projection operator from L'^{B) to the subspace spanned by the 
first N basis functions. Then we find that in Hilbert space notation the A^-th Galerkin 
approximation satisfies 

^ ^-at pNjj^ ^ r e-»it-^)pNKF{U!') ds + e Of . 
Jo 

Here we use to be shorthand for the truncated stochastic convolution 

Of := V A, / e-"(*-^) d(5l Vi . (58) 
i=i 

Hence we obtain for the error of the Galerkin approximation 

Ut-U^ = e-"*(f/o - P^f/o) + J e-'^^'-'^I^KFiUt) - P^KF{UI 
+e (O, - Of) . 

Adding and subtracting the obvious terms yields for the norm the estimate 

Wt-Ul'W < e-"*||f/o-P^?7o|| + IIIP'^/^III r e-"(*-^) II P(f/,) -P(t/f) II rfs 

Jo 

+ \\\K - P^K\\\ re-"(*-^)||P(t/,)|| ds + e \\0t - Of || , 
Jo 
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ds 



where |||P'^i^|||i2 < |||i^|||i2 and sup^gj^ HlP^ii'lllo < oo as a consequence of [3], Lemma 
12.1.4] (cf. the apphcation of this result below). Next, using the Lipschitz and linear growth 
conditions on F, applying Gronwall's inequality, taking the supremum over all t G [0, T] and 
estimating using the bound fl57|) yield 



sup ||f/t-f//^|| < C f ||f/o-P^f/o|| + (l + ||f/o|| + sup 110(1 

t&[o,T] y ^ te[0,T] 

+ C I sup llOt - of II I . (59) 
V*e[o,r] J 

It remains to show that the individual terms in the right hand side converge to zero for 
— )■ oo almost surely. 

• It clearly holds that ||t/o — -P^^o||l2 — and the convergence \\Uo — P'^ Uo\\o — t- holds 
by assumption. 

• Next, as argued above ^1 + \\Uo\\ + sup^gjo.r] is a.s. finite and the compactness 
of the operator K implies jjl-ft' — P^/ir||| -> for — )■ oo, see [3l Lemma 12.1.4]. 

• Finally, the third error term supj^jQ^] jjO^ — Of || vanishes if the Galerkin approxima- 
tions of the Ornstein-Uhlenbeck process O converge almost surely in the spaces 
O([0, T], L'^{B)) and O([0, T], 0(i3)), respectively. This convergence is proven in Lemma 
lAJJ below. 

The proof is completed. □ 

The following lemma contains the convergence of the Galerkin approximation of the 
Ornstein-Uhlenbeck process necessary for proving Theorem 18.11 

Lemma A.l. There exists a sequence > with limTv^oo ^at = such that for all T > 
and all 6 > there exists a random variable Zg with EjZ^p < oo for all p > 1 such that 

sup ||0(-0f|U2 < Zsb]^' 

tG[0,T] 

almost surely. If, in addition, the series Yli^i^f'^i converges in C{B) and the functions Vi 
are Lipschitz continuous with Lipschitz constants Li such that there exists an p & (0, 1) with 
Y^^i ■^f'^'i^ < then it further holds that 



sup ||0(-0f||o < Zsb],-' 

te[o,T] 



almost surely. 



Remark A. 2. Assumptions on the speed of convergence of the series and X^i^i '^'i'^^ 



^ _ ^ *i ^« 

and Xli^i '^i-^'i^ readily yield a rate of convergence for the Galerkin approximation due to 
the definition of the constants in the proof of the lemma. 
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Proof. As in the proof Theorem 18.11 the unspecified norm || • || denotes either the norm in 
LP'iB) or in C{B) and estimates are vahd in both cases. We fix T > 0, p G (0, 1) and a p G N 
with p > 2d /p. Throughout the proof C > denotes a constant that changes from fine to 
fine but depends only on the fixed parameters T,p, p, a and the domain B C M'^. 

Then we obtain for all N,M & N with M < N using the factorization method (cf. 
Sec. 5.3]) similarly to the proof of |8l Lemma 5.6] the estimate 

(Esup,,[o,r] llOf - OfrY^' < C sup,,,,r] (EHF.^^'^r)^/^ , 
where y^^'^^ is the process defined by 

=A/+1 



i=M+l "'O 

In order to estimate the p-th mean of the process y*'^'^ we proceed separately for the two 
cases L^iB) and C{B). 

The case of L'^{B): Due to the orthogonality of the basis functions and employing Holder's 
inequality one obtains 

p/2 

2 



i=M+l 

2(p-2) / X 

i=M+l 



2 




=A/+1 ^ 



t \ p 



Next, as the stochastic integrals in the right hand side are centered Gaussian random vari- 
ables, [H Lemma 5.2^ yields for all t < T 



N \ 2 N , M \P/'i 



\i=M+l / i=M+l ^"^0 ^ 
/ N . „T .p/2 



< 



< 



yi=M+l / j=A/+l 
Af \ P/2 



^2 

Af+1 



^For a centred Gaussian random variable Z it holds EZ^ < p!(EZ^)^'/^ for all p G N. 
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Therefore we obtain for all M,N with M < N 

1/ / ^ \^^^ / °° \^^^ 

(sup,,p,T]Er.^'^1i.) <C[J2 < ^ E ' (60) 

\i=M+l J \i=M+l ) 

where the final upper bound decreases to zero for M — )■ oo by assumption. 

The case of C{B): In this case the estimates get a bit more involved. As p/2 > d/p 
The continuous embedding of the Sobolev-Slobodeckij space W^'^'^iB) into C{B) (cf. [531 
Sec. 2.2.4 and 2.4.4]) and [SI Lemma 5.2] yield the estimates 



suPte[o,r]E||ll ' III; < Csup,g[o,T] /^y^ 1^ _ ^|d+pp/2 ^^^2/ 



Csup,,fo,T] / E|F*^'^(x)rrfx 



B^B k — 1/1 




< CsUPie[0,T] /_ /_ ^^^^ 



/" N p/2 

+ Csup,,fo,T]y^(E|r^''^(x)p) rfx. (61) 

We proceed estimating the two expectation terms in the right hand side. Then we obtain 
for all M < and all x,y E B for the first term 



E|F,^'^(x) 



i=M+l 
N T 

< 5^ AM s-^e-^'^^ds \v,ix) - v,iy)f 

i=M+l 
Af 

< C J2 >^lLf \x - y\''^ (62) 

i=M+l 

for any p G (0, 1) and for the second term 

N t N 

E|r,^^'^(x)|2 < 5^ AM (t-s)-''e-2-(*"'')rfsi;,(a;)2 < >H^h{xf . (63) 

j=A/+l j=M+l 

Next applying the estimates fj63|) and fl62|) to the right hand side of fl6T|) yields, note that 
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pp/2 ~d>0, 



^ 'ill- — <^^<^y-l{j:w) 



i=M+l 



< c([ [ \x - yr^'-'dx dy [ Yl A?A'^)'^'+(sup| A^.(x)^|)'^') 

^•^13 i=M+l ^^'^ i=M+l ' 

, N X 1/2 

sup| ^"^"n 



< c 



i=M+l i=M+l 



for any p G (0, 1). Due to the assumptions of the lemma the two summations in the right 
hand side converge for — )■ oo and thus we obtain for all M, iV G N with M < N the 
estimate 

l/n y oo cx) N 1/2 

(sup,,[o,T]Er,^'^||^) <C sup| Y \'-^(^)'|+ E ^'A'O ' ^^^^ 

^^^'^ i=Af+l i=M+l ' 

where the right hand side decreases to zero for M — oo. 

Overall, we infer from the estimates ( l60l) and ( 164|) that is a Cauchy-sequence in the 
two spaces C([0, T], i^^(-B)) and C([0, T], with respect to convergence in the "p-Xh mean 

and the limit is given by the process O. Moreover, it holds that 

(Esup,g[o,T] llOt-Of f^^'' < C6jv ViVGN, (65) 

where the constant C depends only on p but is independent of and the sequence is 
independent of p and lim7v->oo ^at = 0. As we fixed p G N sufficiently large at the beginning 
of the proof, the result ( 165|) holds for all suffieciently large p G N. Then, however, Jensen's 
inequality implies that ( l65l) holds for all p G [1, oo). Proceeding as in the proof of [HI Lemma 
2.1] using the Chebyshev-Markov inequality and the Borel-Cantelli lemma one obtains that 
there exists for all 5 > a random variable Zs with ElZ^p < oo for all p > 1 such that 

supjg[Q II Ot — of' 1 1 < Zsb]^^ almost surely. (66) 

The proof is completed. □ 
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